Carbon Nanotubes as Cooper Pair Beam Splitters. 
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We report on conductance measurements in carbon nanotube based double quantum dots con- 
nected to two normal electrodes and a central superconducting finger. By operating our devices as 
beam splitters, we provide evidence for Crossed Andreev Reflections tunable in situ. This opens an 
avenue to more sophisticated quantum optics-like experiments with spin entangled electrons. 
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Quantum optics has been an important source of inspi- 
ration for many recent experiments in nanoscale electric 
circuits [1, One of the basic goals is the generation 
of entangled electronic states in solid state systems. Su- 
perconductors have been suggested as a natural source 
of spin entanglement, due to the singlet pairing state of 
Cooper pairs. One important building block required for 
the implementation of entanglement experiments using 
superconductors is a Cooper pair beam splitter which 
should split the singlet state into two different electronic 
orbit als [3, 4]. 

The basic mechanism for converting Cooper pairs into 
quasiparticles is the Andreev reflection in which an orig- 
inally quantum coherent electron pair in the singlet spin 
state is produced at an interface between a supercon- 
ductor and a normal conductor. Conventional Andreev 
Reflections (AR) are local and cannot readily be used 
to create bipartite states [5, 6]. It has been suggested to 
make use of electron-electron interactions 
[ill . [12I . [l3t , sp in filtering [14] or anomalous scattering in 
graphene [15| to promote Cooper pair splitting i.e. the 
Crossed Andreev Reflection (CAR) process. 

In this letter, we show that Coulomb interactions as 
well as size quantization can favor the CAR processes in 
carbon nanotubes. We use a double quantum dot ge- 
ometry where the nanotube is connected to two normal 
electrodes and a central superconducting finger. By op- 
erating our device as a beam splitter (i.e. biasing the cen- 
tral superconducting electrode), we find that there is a 
finite current flowing from the superconducting electrode 
to the left (L) arm and the right (R) arm for a bias voltage 
smaller than the energy gap of the superconductor, which 
demonstrates Cooper pair injection. This subgap current 
is enhanced when we tune the device to the degeneracy 
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FIG. 1: a. SEM image of a typical Cooper pair splitter device 
in false colors with the two biasing schemes sketched. The bar 
is 1/im. A central superconducting electrode is connected to 
two quantum dots engineered in the same single wall carbon 
nanotube (in purple) which bridges between electrodes L and 
R. b. The elementary processes which carry current in the su- 
perconducting (S) state. In addition to the conventional local 
Andreev Reflection process, the Crossed Andreev Reflection 
can occur in which a Cooper pair is split in the two quantum 
dots. The relative probability of each of these processes can 
be inferred from the topology of the beam splitter. 



points of the double dot with the help of capacitively 
coupled side gate electrodes. This enhancement together 
with the dependence of the asymmetry of transport in the 
superconducting state with respect to the normal state 
provide evidence for Crossed Andreev Reflections tun- 
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able in situ^ in contrast to the weakly interacting metallic 
casell6lll7|. These conclusions are supported by theoret- 



ical calculations based on a straightforward modeling of 
the device. 

We use chemical vapor deposition to fabricate our 
SWNTs, which are localized with respect to Au align- 
ment markers by Scanning Electron Microscopy (SEM). 
A SEM picture of a typical device is shown in figure [T^. 
We fabricate our devices using standard e-beam lithogra- 
phy and thin film deposition techniques. We deposit the 
normal and the superconducting contact in one fabrica- 
tion step using shadow evaporation techniques. The nor- 
mal contacts consist of 5nm of titanium followed by 50nm 
of Au or of 70nm of Pd. The central superconducting 
electrode, 80nm — lOOnm wide, is a Al(100nm)/Pd(3nm) 
bilayer. Such a method allows to achieve contact resis- 
tances as low as 30kOhm between normal and supercon- 
ducting reservoir. In addition to the highly doped Si 
substrate with 500nm Si02 which is used as a global 
backgate, we implement two side gates whose voltage 
VsGi and VsG2 can be tuned to control the two different 
parts defined by the central superconducting electrode. 
The spacing between the two normal contacts is between 
600nm and 1.2/im. All the measurements presented in 
this letter have been carried out in a dilution refrigerator 
with a base temperature of SOmK, on one particular sam- 
ple which fulfilled the (stringent) requirements of double 
dot spectroscopy and high enough coupling to the su- 
perconducting electrode. The currents flowing through 
the different arms of the beam splitter are measured via 
the voltage drop across two 2kQ resistors placed in series 
with the device as shown on figure [T^. 

For characterization, we first operate the device as a 
series double quantum dot by setting Vm = and Vs 7^ 

using the bias scheme shown in figure la. Figure 2a 
displays the color scale plot of linear conductance Gl = 
diL/dVs through the right (R) arm of the device as a 
function of side gate voltage 2 Vsg2 and side gate voltage 

1 VsGi- The characteristic 'honeycomb' stability diagram 
of a double quantum dot with rather regularly spaced 
avoided crossings is observed. 

For the remaining of the letter, we operate the device 
as a beam splitter by setting Vm 7^ and l/s = 0. The 
differential conductance corresponding to the left (L) arm 
and to the right (R) arm of the beam splitter have quali- 
tatively the same dependence as a function of Vm, Vsgi 
or VsG2' If VsGi and Vsg2 are tuned out of resonance, 
one typically measures a = diL / dVM shown in figure 
This demonstrates tunneling into a superconductor 
i.e. the current is strongly suppressed for a bias voltage 
below the energy gap A. Such a feature has been used 
very recently to probe the quasiparticle relaxation time 
in SWNTs[18,]. A fitting to the thermally smeared BCS 
density of states gives A = 85/ieV and an electronic tem- 
perature of lOOmi^. Therefore, for Vm < SS/ieF, one can 
only inject Cooper pairs. Note that the energy gap has a 
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FIG. 2: a. Colorscale plot of the differential conductance Gl 
as a function of side gate voltage 1 Vsgi and side gate voltage 
2 VsG2- b. Differential conductance (black squares with lines) 
measured out of resonance in the middle injection scheme at 
l^OmK and OmT. In green square, the differential conduc- 
tance at 44.5mT. In red solid lines, the BCS fit which yields 
an energy gap A of 85/jeV and an electronic temperature of 
lOOmK. 



reduced value with respect to pure Al since it corresponds 
to the minigap of the Al/Pd bilayer [3]. As shown on 
figure [2)3 in green squares, the BCS gap does not appear 
if a magnetic field of 44.5mT is applied perpendicularly 
to the axis of the superconducting finger. This allows 
us to define the normal (N) state of the beamsplitter for 
which we apply a field of 89mT in order to be sure that 
superconductivity of the Al/Pd slab is absent. The su- 
perconducting (S) state is obtained for OmT. 

The colorscale plot of figure 3a displays Gl for Vm = 
AQjieV as a function of Vsgi and Vsg2 for a specific anti- 
crossing (ACl) in the S state. In contrast to the off 
resonance case shown in figure 2b, a relatively high sub- 
gap current Is ^ 0.1 /at can flow into or from the su- 
perconductor even though Vm is smaller than the energy 
gap at ACl. The presence of a subgap current can only 
be understood if Andreev reflections are taking place at 
resonance. In order to characterize the type of Andreev 
process occurring near the anticrossing, we show in figure 
3b the variations of Gl = diL/dVM and Gr = diR/dVM 
measured simultaneously along the yellow arrow of the 
colorscale plot in figure 3a. In the N state (black solid 
lines), Gl and Gr display two peaks corresponding to 
the bonding/antibonding states of the double dot. Due 
to an asymmetry of the coupling of each dot to the L(R) 
reservoirs, the height of each doublet of peaks is differ- 
ent. In the S state represented in red solid lines, the 
peak height is reduced, but there is still a finite con- 
ductance for the L(R) arm which follows essentially the 
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resonances observed in the N state. A similar feature 
is shown in figure 3c for another anti-crossing (AC2) for 
which the asymmetry between Gl and Gr is shghtly dif- 
ferent. Note that the curves in the N state have been 
scaled down by 1/3 for the sake of clarity. 
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FIG. 3: a. Focus on a specific region for the differential con- 
ductance as a function of Vsgi and Vsg2 (ACl). The 
yellow arrow indicates the direction in which the linear scan 
of panel b. has been taken, b. Measurements of the con- 
ductance in the normal state (black solid lines) and in the 
superconducting state (red lines) for Or and Gl along the 
direction of the yellow arrow. For the sake of clarity, the N 
state conductances have been multiplied by 1/3. The model 
calculations are in green solid lines for the S state and in or- 
ange solid lines for the N state. In blue dashed lines, the CAR 
probability, c. Similar graph as in b. for anticrossing 2 (AC2) 

The elementary processes which contribute to the sub- 
gap current in our double quantum dot setup at the de- 
generacy point are schematically sketched in figure lb. 
The tunnel rates to the S,L and R contacts are respec- 
tively T^, Tl and Tr. The initial state in the supercon- 
ductor S is represented in blue and the final state in the 
normal metals N, in red. In process 1, two particles are 
transferred to the same reservoir and the corresponding 
probability is proportional to r| ^. In processes 2 and 3 
which are equivalent, one particle is transferred to each 
reservoir. The corresponding probability is therefore pro- 
portional to VlTr. While the beam splitter geometry 
imposes the above general form for the Andreev tunnel- 
ing probabilities, our theory (which we discuss in detail 
in the EPAPS) allows to give an absolute value of the 



proportionality constants of the AR and the CAR, by 
accounting for higher processes. 

From the anti-crossing of figure 3a and the full stabil- 
ity diagram of figure 2a, we determine all the important 
parameters of the double dot using non-linear transport 
similarly to what is done in ref. 20 for example. The 
corresponding plots obtained from our theory at T = 
are given in green solid lines for the S state and in or- 
ange solid line for the N state for two anticrossings ACl 
and AC2. Each anticrossing is characterized by the set of 
parameters {U l,Ur,V 12,^ l,^ r,^ sl,^ sr}, Ul,r being 
the onsite charging energy, the coupling between the 
two dots, Tl^ Tr the coupling of each dot to its reservoir 
and TsL, ^sr the coupling of the superconductor to each 
dot. Although not all of these parameters can be deter- 
mined independently, the fitting establishes the following 
hierarchy of energy scales: Ul,r ~ ImeV, ^ lOO/ieT^, 
Tl^Tr ^ 100/ieV and Tsl^^sr ^ 10/ieV. For both 
arms and both anticrossings, the calculated current in 
the superconducting state is in good agreement with our 
experimental findings, which implies an important con- 
tribution from CAR processes, as shown by the dashed 
lines in figure 3b and 3c. For ACl, for example, the fit 
allows us to extract a contribution of split Cooper pairs 
up to 35% for Gl and 55% for Gr. 
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FIG. 4: Top panel : Efficiency of the Cooper pair splitting for 
the device. The solid green line corresponds to the asymmetry 
variation predicted by the interacting theory. Bottom panel 
: Asymmetry in the superconducting state versus asymmetry 
in the normal state for 31 out of the 35 anti-crossings studied 
(red closed squares) and 11 resonance peaks away from the 
ant i- crossings (black stars). The dashed curve corresponds to 
o^s = c^%. The solid dashed line corresponds to as = cxn- 
The green line is the full theory (see EPAPS) and the pink 
line is the full theory without the CAR process. 

Now, we exploit the asymmetry of transport between 
the L and R arms to provide further evidence for the 
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Cooper pair splitting. In the S state, it is defined as 
^5- = Gl/Gr at OmT, and in the N state, as = 
Gl/Gr at 89mT, both taken at resonance. One can al- 
ready anticipate from the elementary processes of figure 
lb that the presence or the absence of processes 2 and 3 
should be observable directly from the dependence of 
with respect to a at. At the degeneracy point, both the 
local and non-local Andreev processes acquire a Breit- 
Wigner-like form at resonance (see EPAPS). In the limit 
>> Tl^r >> r^, obtain for the conductance GL(^Ry. 
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where is the tunnel rate to the superconductor renor- 
malized by Coulomb interactions. The presence of the 
crossed term F^^Fi^ in equation (1) implies that: 
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Away from the degeneracy points, when one of the dots 
is blocked, only local AR are present and the term pro- 
portional to F^Fi^ in equation (1) is absent. In this case, 
as is just the square of a at. The bottom panel of figure 
4 displays as versus aN for all the anti-crossings which 
we have measured, in red squares, and resonances away 
from the anticrossings, in black stars. The 31 red squares 
correspond to 7 different anticrossings. The error bars 
correspond to the systematic error made when determin- 
ing as. The black stars are obtained when one of the 
two dots is Coulomb blockaded (single resonance case). 
Within the error bars, they fall onto the dashed curve 
which corresponds to the universal parabola as = aj^ as 
expected. The red points fall onto the universal dashed 
line as = a^ [21]. The observed contrast between the 
quadratic behavior followed by the black stars and the 
linear behaviour followed by the red squares proves the 
Cooper pair splitting action of our device. Note that the 
full theory (i.e. without the approximations leading to 
Eq. (1)) predicts a dependence of as versus aN which is 
marginally different from these behaviors in our parame- 
ter range as shown by the green and pink curves in figure 
m The universal linear behavior arising from the Cooper 
pair splitting is therefore a very useful diagnosis tool for 
a wide range of parameters of our beam splitter. The 
efficiency of the beam splitter can be defined as the ratio 
2Tcar/{ARl + ARr + 2Tcar) - 2/(2 + 1/as + as) of 
the split current to the total current. It is displayed on 
the top panel of figure [H For a a at ranging from 0.5 to 
1, the efficiency of Cooper pair splitting is close to 50%, 
showing that almost half of the Cooper pairs ffowing out 
of the superconducting finger are split at the SWNT. 

In conclusion, we have shown that carbon nanotube 
double quantum dots can be used as tunable Cooper pair 
beam splitters. The specific advantage of beam splitters 
based on double quantum dots is the possibility of further 



processing the electron states e.g. by spin filtering for 
EPR type experiments. 

During the course of writing this paper, we became 
aware of a similar study using InAs nanowires by Hof- 
stetter et al.|22|]. 
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